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Abstract 

We estimate CP violation for several experimental facilities study- 
ing neutrino oscillations. We also estimate the probability of ia, to u e 
conversion, using values of the parameter #13 within the known lim- 
its, in order to suggest new experiments to measure CP violation for 
neutrinos moving in matter. 



1 Introduction 

The study of CP violation (CPV) is essential for understanding weak inter- 
actions. Almost half a century ago CP violation in weak interactions was 
found in the decay of K\ into 7r + + 7r~ pQ and 27r° [2J, with branching ratios 
of the order of .001. The decay K\ — > 7r° + v + v is almost entirely CP 
violating [3J but requires accurate determination of the CKM matrix jl] and 
accurate measurements. See Ref[5] for a review of this experiment and ref- 
erences. There have many other studies of CP asymmetries in weak decays: 
see Ref jH] for a recent study of B radiative decay with references to earlier 
work on CP violation in various weak decays. 
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In our present work we study possible CP and CPV measurements using 
neutrino oscillations. In recent years there have been a number of experimen- 
tal studies of neutrino oscillations using neutrino beams from accelerators and 
reactors, and a most important objective of these experiments is the mea- 
surement of CP violation (CPV). In our present work on estimating CPV we 
use parameters for the baseline and energy corresponding to MiniBooNE[7], 
JHF-Kamioka [8], MINOS [9], and CHOOZpJJ], which are on-going projects. 

There have been many recent studies of CP and T symmetries via neu- 
trino oscillations for future facilities, e.g., see RefsfTTl [12], which also give 
references to earlier publications, and the ISS report [13J on future neutrino 
facilities. One possible future facility for studying CPV and the 5c p pa- 
rameter is the LBNE Project, where neutrino beams produced at Fermilab 
would have a baseline of L ~ 1200 km, being detected with deep underground 
detectors [T4"l IT5] . With the methods used in the present work, described be- 
low, predictions of CPV with the baseline and energies of the LBNE Project 
have recently been made for 8c p from 90 to degrees [T6]. The angle 6*13 is 
not well known, and will be measured by the Daya Bay experiment in China. 
First, we use Sep = 90° and two values of # i3 to explore the dependence of 
CP and CPV neutrino oscillations on this parameter for all on-going projects. 
We then calculate CPV for JFK-Kamioka baseline with E=.48 GeV, which 
has a large CPV for sin9i3—0.19 (Sec. 3), to estimate CPV for values of #13 
expected to be found (Sec. 4). 

A major complication for the determination of T, CP, and CPT violation 
is the interaction of neutrinos with matter as they travel along the baseline. 
These matter effects have been studied by a number of theorists. See, e.g., 
Refs[T71 Q21 US]- The main objective of the present research is to estimate 
matter effects for CPV. for the MiniBooNE, JHF-Kamioka, MINOS, and 
CHOOZ facilities. 

For the basic interactions, which are CPT invariant for local theories, 
CP violation also implies T violation. Our present research is an extension 
of our recent work on T reversal violation[2U]. In that study we used the 
formalism of Ref[2T] for v e <H- TRV, and that of Ref[22] for v e — > 
conversion probability to calculate the effects of neutrinos moving through 
matter. In the present work we use the notation and formalism of Jacobson 
and Ohlsson[23], who studied possible matter effects for CPT violation. 

CP violation in the a — b sector is given by the transition probability, 
denoted by V[y a — > z/5), for a neutrino of flavor a to convert to a neutrino 
of flavor b; and similarly for antineutrinos v a -,^b- The CPV probability dif- 
ferences (note that the C operator changes a particle to its antiparticle) are 
defined as 

APf fe P = V{v a ^v h )-V(v a ^v h ) . (1) 
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In our present work we study V{y^ — > z/ e ) and V(u^ — >■ P e ), since the 
neutrino beams at MiniBooNE, JHF-Kamioka, and MINOS, as well as most 
other experimental facilities, are muon or anti-muon neutrinos. 



2 Transition Probability Viv^ — > v e ) 

In this section we review the derivation of the probability of a muon neutrino 
to convert to an electron neutrino, Viy^ — > z/ e ), using the notation of Ref[23j. 
We then make an estimate of the transition probalilities for sample acceler- 
ator and reactor experiments. Although at the present time no experiments 
for CPV are possible, this can serve as a basis for future experiments. In the 
next section we give somewhat more accurate calculations for CPV for the 
same set of experimental facilities. 

As in Refs[2"Tj 123"] we use the time evolution matrix, S(t,t ) to derive 
the transition probabilities. For neutrino oscillations the initial neutrino 
beam is emitted at time to, usually taken as to = 0, and the neutrino or 
converted neutrino is detected at baseline length = L at time=t. Since the 
neutrinos move with a velocity near that of the speed of light, at the end of 
our derivation we take t — to — > L, with the units c=l. 

Given the Hamiltonian, H(t), for neutrinos, the neutrino state at time = 
t is obtained from the state at time = to from the matrix, S(t,t Q ), by 

\v{t)> = S(t,t )\u(t ) > (2) 
^S(Mo) = H(t)S{t,to) . (3) 

Neutrinos (and antineutrinos) are produced as u a , where a is the flavor, 
a = e, fi, t. However, neutrinos of definite masses are u a , with a = 1,2,3. 
The two forms are connected by a 3 by 3 unitary transformation matrix, U : 
v a = U u ai where z/ a , v a are 3x1 column vctors and U is given by (sm% = s^-) 



/ C12C13 S12C13 si 3 e '" GfJ 

U = -S12C23 - Ci2S23Sl 3 e 4<5cp C12C23 - s 12 s 2 3Si 3 e tScp S23C13 
\ S12S23 — Ci2C2zS\ 3 e %&CP — C 12 S 2 3 — Si2C 23 Si3e j5cp C23C13 

similar to the CKM matrix for quarks. We use the best fit value [TT] S23 = 
0.707. #13 is not well known. We use S12 = 0.56. We use S13 = 0.19, consistent 
with a recent analysis[l2], and Si3=0.095, to determine the dependence of CP 
and CPV on this parameter, which is not well known. The CP phase 5c p is 
also not well known. For simplicity we choose 5qp = vr/2, and calculate the 
dependence of Viy^ — > z/ e ) on 5qp, as discussed below. 
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In the vacuum the S(t, to) is obtained from 



Sab{t,t ) = J2 U ^ eX P iEj(t ~ t0)U bj ■ ( 4 ) 
3=1 

Since neutrino beams in neutrino oscillation experiments travel through 
matter, and the main neutrino-matter is scattering from electrons, we must 
include potential, V = \f2Gpn e „ for neutrino electron scattering in the earth: 
where Gf is the universal weak interaction Fermi constant, and n e is the 
density of electrons in matter. Using the matter density p=3 gm/cc, the 
neutrino-matter potential is V = 1.13 x 10 -13 eV. 

The transition probability V{v^ —> u e ) is obtained from Si 2 , with 



7>K -> ve) = \Si 2 \ 2 = Re[S 12 } 2 + Im[S 12 ]\ with 

S12 = c 23 P-is 23 ae-* Scp A (5) 

a = s 13 (A-sj 2 5) (6) 

5 = 5m\ 2 /{2E) (7) 

A = 5m\J(2E) (8) 

A ~ f(t)I a * (9) 

I a * = [ dt' a* it') fit') (10) 
Jo 

a(t) = cosut — icos26sinut (11) 

/(*) = e" 4At (12) 

2co = ^5 2 + V 2 - 25Vcos{26 l2 ) (13) 

(3 = —isin29sinujL (14) 

A = A-{V + 5)/2 (15) 

sin29 = S12C12— , (16) 



where the neutrino mass differences are 8m\ 2 = 7.6 x 10 5 (eV) 2 and 8m\^ = 
2.4 x 10~ 3 (e1/) 2 . Note that 5 < A, and t L for u„ ~ c. From Eqs. ffTTTfHl) : 

Re[Si 2 ] = s 23 a[cos((A + 8cp)L)Im[I a *\ - sin((A + 5cp)-^)-Re[/ Q *] 
Jm[S'i2] = —c 23 sin29sinuL — s 2 3a[cos((A + <5cp)L)i?e[/ aJk ] 

+sm((A + 5 CP )L)/m[/ Q ,]] (17) 
Viy^^rv^) ~ (c23Si2Ci 2 (5/ w ) s ^ w ^) 2 + (^sSissmAL) 2 ( 18 ) 

+ (S23S1 3 (C0SAL - l.)) 2 + 2S13S12C12S23C23 

(S /u)sinuL(cosAL — 1) . 
From Eq (jl8p we obtain the results for Viv^ — > v e ) shown in Fig.l. 
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Figure 1: The ordinate is Viv^ — > v e ) for MINOS (L=735 km), 

MiniBooNE(L=500m) , JHF-Kamioka(L=295 km), and CH00Z(L=1 km). 
Energy= E in GeV. Solid curve for s 13 =.19 and dashed curve for si 3 =.095 
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Although these results are only for muon to electron neutrino conversion, 
they can provide guidance for future experiments on CPV via <H- v e os- 
cillation. Note that in Ref[8] Viy^ — > v e ) was calculated for the 295 km 
JHF-Kamioka project for E=0-2 GeV, and our calculation based on the the- 
ory developed in Refs. [23ll2~T] . finds that with similar parameters Viv^ — > u e ) 
is in agreement for E=. 4-1.0 GeV with this earlier estimate. Since there is 
uncertainty in the value of 5cp, we calculated Viy^ — > v e ) for Sep = 0. We 
do not show the results, as they are almost the same as shown in Fig. 1. 



3 CP Violation AV^f 

In this section we shall extend the derivation of the transition probability 
Viy^ — > z/ e ) of the previous section to derive the CPV probability 

= |5 12 | 2 -|5 12 | 2 (19) 

with 5*i2 defined in Eq(8) and 

3i2 = c 23 (3 - is 23 ae lScp A , (20) 

with f3 = (3(V ->■ -V) and A = A(V = -> -V). For example (see Eqs(16,18)) 
2uj = ^5 2 + V 2 + 2SVcos(29 12 ) and A = A + (V - 5) /2. Using conservation 
of probabiltiy|23], \A\ 2 = \A\ 2 . With 5 CP = tt/2, e^l + ^ lScp = 

W C J = C 2 3(|/3| 2 -|^| 2 )-2c 23S2 3a(/mH/3A*]-/m^]). (21) 

From Eq fl2T|) . the definitions in the previous section, defining s = sin(ujL), 
c = cos(ujL), and using 5c p = vr/2 one finds 

2 —2 

A7f e P = c 2 23 s 2 12 c 2 12 5 2 (^-^) + 2c 23 s 23 s 12 c 12 s 13 5(A-5s 2 12 ) (22) 
r uj z tu z 

,s, -r r . A — ucos29 s, = T .. A — ujcos29 . 

- c - cosAL — — + - c - cosALW — . 

uj A 2 — uj 2 UJ A 2 — uj 2 

The results for AV^f are shown in Fig. 2. Note that the largest values for 
CPV are for CHOOZ, with a small baseline and low energy. With Kamioka 
parameters, CPV would also be a few percent, which might be measurable. 
However, for experimental tests of CPV one needs both and beams 
with the same parameters. Perhaps this will be possible in the future. 

As has been stated in many publications, in vacuum AV^'f is given by 
APjg, and both vanish if 5qp = 0. However, with matter effects (V ^ 
0), AV^f and APj e , must be treated separately, as we have done. The 
magnitude of 5c p is important for predictions, and it is expected that in the 
future it will be determined with greater accuracy. 
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Figure 2: The ordinate is AP(z/ M ->■ i/ e ) for MINOS (L=735 km), 
MiniBooNE(L=500m) , JHF-Kamioka(L=295 km), and CH00Z(L=1 km). 
Energy=E in GeV. Solid curve for s 13 =.19 and dashed curve for s 13 =.095 
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4 AV^p For JHF-Kamioka with E=0.48 GeV 
and sin6i3=0.0 to 0.19 

From Fig. 2 one sees that with the value sin6\ 3 = 0.19 CPV reaches a value 
of over 4%, for the JFK-Kamioka project, which could be detected if beams 
of both neutrinos and antineutrinos were available. Since the value of sinOiz 
is not known at the present time, we use the energy and baseline values of 
E=0.48 GeV, L=295 km with sin0 13 from 0.0 to 0.19[12]. Using Eq.(f22]), 
with sinOiz a variable, one obtains the results are shown in Fig. 3. 
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Figure 3: AP(z/^ ->- v e ) for JHF-Kamioka(L=295 km), E=0.48 GeV, as a 
function of sindn 

Note that since the second term in Eq. (l22|) is dominant, the dependence 
of AP(f /1 — > v e ) on sin8i3 is almost linear. 
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5 Conclusions 



We have estimated CP violation for a variety of experimental neutrino beam 
facilities. No experiments are possible now to test CPV via neutrino oscil- 
lations, since beams of both neutrino and antineutrino with the same flavor 
would be needed, with parameters chosen for a CPV of 1% or more to make 
the experimental measurement possible. Our results should help in planning 
future experiments. Note, however, that our results depend on the value of 
S13 and Sep, which are not well known, and we have used two values for s 13 
to determine the dependence of CP and CPV on this perameter (see Fig. 2), 
and estimated CPV as a function of S13 for the JFK-Kanioka baseline of 295 
km and energy E=0.48 Gev, with the CPV probability over 4% for s 13 =0.19 
(see Fig. 3). 
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